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ABSTRACT 

Derivation  of  two-dimensional  equations  for  elastic  deformations 
of  laminated  anisotropic  plates,  based  on  the  assumption  of  a Kirchhoff- 
distribution  of  primary  strains,  in  conjunction  with  the  use  of  the  Castigliano 
variational  principle  for  stresses,  or  a variational  principle  for  stresses 
and  displacements.  A discussion  is  given  of  the  relation  of  the  present  work 
to  some  earlier  work. 


Note  on  the  Effect  of  Transverse  Shear  Deformation 

t 

in  Laminated  Anisotropic  Plates 
By  E.  Reissner 

Introduction.  The  recent  work  by  Argyris  et  al.  [l]  and  by  G.  A. 

Cohen  (2]  on  the  approximate  determination  of  transverse  shear  stress 
coefficients  for  laminated  anisotropic  plates  (and  shells)  has  shown  the 
relative  complexity  of  the  task  of  extending  the  earlier  analysis  of  the  effect 
of  transverse  shear  in  isotropic  homogeneous  [3]  as  well  as  sandwich-type 
plates  [4]  to  the  case  of  laminated  anisotropic  plates.  The  point  of  departure 
of  the  following  considerations,  which  were  motivated  by  Cohen's  paper,  is 
a brief  note  [6]  on  an  earlier  rational  approach  to  the  subject  of  transverse 
shear  stiffness  in  laminated  anisotropic  plates,  which,  for  brevity's  sake,  was 
limited  to  the  case  of  plates  with  such  symmetry  as  to  allow  a separate  treat- 
ment of  transverse  bending  and  stretching. 

A study  of  the  contents  of  [2]  made  it  appear  that  the  results  in  this 
paper,  although  quite  different  in  appearance  and  derivation,  should,  for 
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the  symmetrical  case,  be  in  fact  equivalent  to  the  results  reported  in  [6], 

This  being  the  case,  the  present  brief  paper;  although  not  presenting  numerical 
results  or  procedures,  is  thought  to  be  of  some  interest  from  the  point  of  view 
of  the  development  of  computer  methods  in  applied  mechanics. 

Briefly,  the  present  work  as  well  as  the  work  in  (2)  and  (hj  starts  out 
with  a distribution  of  primary  stresses  in  laminated  plates  corresponding  to 
a linear  (Kirchhoff)  distribution  of  primary  strains.  The  second  step  in  the 
three  papers  is  a determination  of  the  associated  statically  consistent  trans- 
verse shearing  stresses.  The  remaining  step  is  the  use  of  an  energy 
principle  for  the  utilization  of  these  expressions  for  primary  and  transverse 
stresses  for  the  purpose  of  obtaining  two-dimensional  constitutive  equations 
involving  transverse  shear  deformations  and  transverse  shear  stress 
resultants.  In  [6]  this  third  step  is  executed,  not  quite  as  simply  as  it  might 
have  been  done,  through  an  application  of  the  writer's  variational  principle 
for  stresses  and  displacements  [5].  In  Cohen's  work  f2]  this  third  step  is 
executed  through  an  application  of  the  principle  of  minimum  complementary 
energy  (variational  principle  for  stresses).  In  the  present  reconsideration 
of  the  problem  two  types  of  results  are  presented.  These  are  associated, 
respectively,  with  the  variational  principle  for  stresses,  or  with  the  variation.-.l 
principle  for  stresses  and  displacements.  It  is  found  that,  while  Cohen's  method 
[2]  depends  on  use  of  the  stress  principle,  his  results  appear  to  be  equivalent 
to  the  results  in  f6]  which  depend  on  the  use  of  the  stress  and  displacement 
principle. 
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Stress  Strain  Relations  and  Expressions  for  Stresses.  We  consider 
a plate  with  midplane  coordinates  x , x2>  and  with  thickness  coordinate  z. 

We  designate  midplane  parallel  stresses  by  o , O^,  (with  representing 
the  shear  stress  oj2  = ).  and  we  designate  transverse  shearing  stresses 

by  r , T . We  assume  that  the  deformatlonal  effects  of  transverse  normal 
stress  - which  would  here  be  designated  by  T ^ - are  negligible,  and  we  assume 
that  the  relevant  stress  strain  relations  for  the  laminated  anisotropic  plate  are 
of  the  form 


o.  = E..e. 
i *1  J 


TX  = GxuRu 


(la.b) 


with  the  coefficients  E^  and  being  known  functions  of  z. 

Given  the  stress  strain  relations  (la),  we  now  introduce  as  a basic 

* 

assumption  for  the  analysis  to  be  undertaken  that  the  distribution  of  face- 
parallel  stresses  may  be  approximated  in  the  form 

■ Vi  *V  ■ 

with  Xj  and  (,  being  independent  of  z. 

Equations  (2)  imply  as  relations  for  stress  couples  and  stress  resultants. 


“•  ■ Vj  ♦Vj  • N>  ’ cu"j + b<jcj  • 


where 


(Bj.,  Cj.,  D^)  = /_cCl.  z.  zZ)Eijdz  . 

We  agree  to  write  the  inverted  form  of  (3)  as 

-I  ■ Dij1vc<VNi  • ‘j  ■c«MjtB«Hj  • 


(3a.  b) 


(4) 


(5a,  b) 


Introduction  of  15)  into  equation  (2)  leads  to  approximate  expressions  for  face- 
parallel  stresses  in  terms  of  stress  couples  and  stress  resultants  of  the  form 


* 9 0 6 


o 
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"l  ■ "uwj  • 

where 

h.j  * <*  * SjXk  • r„  = (zckj  * %'>E,k  • 


<*>) 


(7) 


Having  equations  (6),  we  uae  as  associated  approximate  expressions  for 

the  transverse  stresses  the  equilibrium  expressions 
z z 

T1  ■ -•r-c<',l.lt,’).2|d*  • T2  ■ -0,,J.14°2.Z»‘U  • <»> 

in  conjunction  with  the  conditions  of  vanishing  T ^ for  z = +c.  It  is  evident  from 

(6)  and  (8)  that,  with  suitable  coefficients  K and  J,  these  expressions  can  be  written  as 

T.  ~ K\  + J.  tllN,  . (9) 

X XjU  j.u  Xju  j.u 

Having  equation  (9)  we  now  introduce,  in  addition  to  the  resultants  N and 

j 

the  couples  M^,  transverse  shear  resultants  Q^,  and  we  observe  that  the  set 
Nj.  subject  to  the  five  two-dimensional  equilibrium  equations 

Ml,l+M3,2  = Q1  * M3,  lfM2,2  = °2  * Ql,l+Q2,2  =0  • (10a) 

N1,I*N3,2  3 0 ' N3,l+N2,2  3 0 * (,0b> 

We  may  use  equations  (JO)  in  order  to  express  the  T ^ - instead  of  in  terms 
of  the  twelve  y and  ^ - in  terms  of  ten  quantities  only,  these  being  the 
two  transverse  resultants  Q^,  and  two  sets  of  four  quantities  Rj.  and  Sj,  defined  by 
R,  = M,  , - Mj  2 , R,  = M j 2 , R,  * M,  , - M2  2 , R4  = Mj  ,,(11*1 

S1  ’ N1.1-NJ,2  • S2  ’ N,,2  • S,  ■ N).l-N2.2  • S4  * N2. 1 • 


The  resulting  expressions  for  may  be  written  in  the  form 

TX  = AXuQu  + l\L  R£  + LXESE  ’ 


(12) 


with  the  problem  of  the  appropriate  treatment  of  the  quantities  R and  S being 
of  particular  interest  in  what  follows. 
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Bending  and  Stretching  Uncoupled.  It  is  sufficient  for  present  purposes 
to  limit  further  discussion  to  the  uncoupled  problem  of  transverse  bending,  which 
is  given  when  all  the  coefficients  and  in  (1)  are  even  functions  of  the 
thickness  coordinate  s.  We  obtain  the  equations  for  uncoupled  transverse 
bending  by  setting 

«j  = 0 ' CiJ  ’ ° • Fij  = ° • JXjM  = 0 ■ SE  " 0 * <13> 

with  the  remaining  plate  problem  now  consisting  of  the  standard  equilibrium 
equations  (10a),  in  conjunction  with  approximate  two-dimensional  constitutive 
equations  and  strain  displacement  relations,  with  these  having  to  be  consistent 
with  the  remaining  portions  of  equations  (Z),  (3)  and  (12). 

Use  of  the  Principle  of  Minimum  Complementary  Energy.  We  have  as 
three-dimensional  strain  energy  density  expression  in  terms  of  stresses,  consis- 
tent with  die  stress  strain  relations  (la,b) 

w = . (ui 

Leaving  aside  boundary  condition  considerations,  we  may  obtain  a 
system  of  approximate  two-dimensional  constitutive  equations  by  utilising  the 
variational  equation 


V 

ft  JJ/  _cWdsdx1dx2  = 0 , 


in  conjunction  with  the  equilibrium  stress  expressions 


* HuMj  • 'x  ■ KXjuMj.u  ■ 


( 16a,  b) 


subject  only  to  the  supplementary  two-dimensional  equilibrium  equation 


Ml,  11  + 2M3,  12  + M2.22  = 0 * 


r 


,iXXuQxQu  ,ou. 

+ (M1,  1 + M3,2  ' Q\),pl  + (M3,  ] + M2.2  ' Q2)<<52Jdxidx2  = 0 * (22) 
The  resulting  equations  are 


D..  M.  - <p, 
ij  J •. 


d:!m.  = f/»-  , 

2j  j 2,2 


DM.  - + (p_ 

3j  J 1.2  v2 , 1 


(23) 


together  with 


X Q = w + to 
1U  U .1  1 


X2MQU  = W,2+,P2 


(24) 


Plates  for  which  some  of  the  ^ and  Z^. ^ are  non-vanishing. 

In  considering  the  variational  equation  (15)  for  this  general  case,  we  must  take 
account  of  all  or  some  of  the  terms  with  Rj,  in  equation  (20),  in  a consistent 
manner.  This  may  be  done  without  additional  Lagrange  multipliers  by  writing 
the  terms  R^,  in  (20)  in  terms  of  derivatives  of  the  M.  in  accordance  with 
equations  (11a)  and  by  then  using  the  one  constraint  equation  (17).  Alternately, 
we  may  introduce  equations  (11a)  as  constraint  equations,  through  additional 
Lagrange  multipliers  tj, . If  we  adopt  the  latter  procedure,  we  obtain  as  a 
generalization  of  the  system  (23)  and  (24)  the  following  set  of  constitutive  equations 


-1 


D,  . M.  = <P,  . + ♦,  , + , 

lj  J 1,1  1,1  2,2 


2j  j 2,2  3,2  4,1 


(25) 


D'M.  = <P.  , +<P,  . - ♦.  , + , 

3j  j 1,2  2,  1 1,2  3,1 


together  with  the  two  relations 


x*pQp  + Y^oRn  = W,r.  + «V 


(26) 


and  the  four  relations 
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YvnQv  + zroRE  3 *n  * (27) 

We  note  that  equations  (25)  to  (27).  together  with  the  equilibrium 
equations  (10a)  and  the  defining  relations  (11a),  are  a system  of  sixteen 
equations  for  the  nine  stress  measures  Q^.  R£  and  seven  displacement 

measures  w,  tj--  W w«  eliminate  the  quantitites  Rj,  and  we  will  be  left 
with  altogether  eight  equations  for  the  M..,  Q^,  w and^.  We  refrain  from  listing 
these  equations  and  limit  ourselves  to  noting  the  important  fact  that  the  presence 
of  the  geometrical  variables  in  (27)  means  that  it  will  now  not  be  possible  to 
take  account  of  the  effect  of  transverse  shear  deformations  by  rrwans  of  constitu- 
tive equations  of  a form  as  simple  as  the  set  (24). 

A supplementary  approximation.  It  is  possible,  formally,  to  use  the 
variational  equation  (15)  for  the  derivation  of  approximate  constitutive  equations 
without  requiring  that  the  transverse  shear  stress  distributions  in  (18)  are  in 
exact  three-dimensional  equilibrium  with  the  stresses  0^  in  (16a)  although,  in 
this  case,  one  then  no  longer  uses  the  principle  of  minimum  complementary 
energy  but  rather  an  ad  hoc  modification  of  this  principle.  One  such  non- rational 
modification  is  to  ignore  the  fact  that  the  quantities  R^.  in  (18)  are  defined  in 
terms  of  the  Mj  . in  accordance  with  (11a)  and  to  treat  the  R^  as  supplementary 
free  parameter  functions.  When  this  is  done  there  is  then  no  occasion  for  the 
introduction  of  the  Lagrange  multiplier  functions  and  therewith  the  constitu- 
tive equation  system  (25)  to  (27)  becomes  simplified  by  the  disappearance  of  all 
♦£ -terms.  This  fact  is  of  particular  significance  insofar  as  equations  (26)  and 
(27)  are  concerned,  because  it  is  then  possible  to  solve  (27)  in  the  form  R = 

‘zravvaQv  and  lo  wrtte  ln  the  want«d (xyp  - v^z^y^kj  . 
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} 


It  appears  that  the  results  reported  in  {Z],  when  specialized  to  the  case 
of  uncoupled  stretching  and  bending,  are  in  fact  equivalent  to  what  is  obtained 


in  this  particular  fashion. 

Use  of  a Variattonal  Principle^ 
the  variational  equation 

1°1  4 “z,*"*  4 <"1,2  * u2. 1>°3  4 •«!,.  4 1»T, 

+ (u  + u )T  - W)dz  dx  dx  = 0 , 

c,  p Ti  Z,  6 b 1 c 


(Z8) 


in  place  of  equation  (15),  is  known  [5]  to  be  consistent,  without  a priori  constraints 
on  stress  distributions  and  with  or  without  independent  assumptions  concerning 
displacements.  Consequently,  if  "we  depart  from  (28),  we  may  now  use  the  strain 
energy  density  expression 
-1, 


(29) 


= «Dl'j  MlMj  4 \oQ"Q° 4 2 ykoqi.ro  4 ZEORt  Rn>  • 

without  considering  (11a)  as  constraint  equations  for  the  . Furthermore,  if 
we  introduce  as  approximating  expressions  for  displacements,  consistent  with 

the  stress  expressions  <7,  = E..zx., 

i ij  j 

UX  = **X(*rV  ' u*  = w(x|»x2)  ' (30a,  b) 

we  then  have  in  equation  (28) 
c 

JcS.lV**'  +(U2,z  + UZ,2)T2ldz  = Vl.  1M1  +<P2.2M2 

+ <*1 .2  +V,2,1)M3  +(<P1  +w  j)0,  2)Q2  * (31> 

and  the  Euler  equations  of  the  remaining  variational  problem  come  out  to  be  the 
three  two-dimensional  equilibrium  equations  (10a),  together  with  constitutive 
equation  of  the  form  (23),  in  conjunction  with  an  abbreviated  system  of  trans- 
verse shear  equations 
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